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^ : Abstract 



We solve the Chapman-Kolmogorov equation and study the exact splitting probabilities of the 
general stochastic process which describes polymer translocation through membrane pores within 
the broad class of Markov chains. Transition probabilities which satisfy a specific balance constraint 



-C provide a refinement of the Chuang-Kantor-Kardar relaxation picture of translocation, allowing us 

to investigate finite size effects in the evaluation of dynamical scaling exponents. We find that (i) 

previous Langevin simulation results can be recovered only if corrections to the polymer mobility 
cd; 

exponent are taken into account and that (ii) the dynamical scaling exponents have a slow approach 

to their predicted asymptotic values as the polymer's length increases. We also address, along with 

S : 

strong support from additional numerical simulations, a critical discussion which points in a clear 
way the viability of the Markov chain approach put forward in this work. 
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I. INTRODUCTION 

The phenomenon of polymer translocation through membrane pores has received a great 
deal of attention in recent years [lj. Important issues are related to the translocation of 
complex biomolecules in the methabolism of living cells [2|, and, on the technological fore- 
front, to developments in the fields of targeted drug/gene delivery [3( and DNA sequencing 

It has been known, however, that there is not a universal mechanism for polymer translo- 
cation. A case-by-case analysis is necessary to conclude if translocation is related to specific 
features of the biochemical environment which surrounds the membrane or to the existence 



of biomolecular motors, as the ones found in mitochondria 



a. 



Taking a look at the volumi- 



nous literature on the subject, one finds studies of polymer translocation driven by chemical 



14. 



or electric potential gradients 
simply unbiased translocation 
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12] , chaperone-assisted rectified brownian motion [13( or 



14 



15| 



The case of unbiased translocation, where a polymer is let to diffuse through a mem- 
brane pore as the sole consequence of thermal fluctuations, is the ideal starting point for 
the investigation of more sophisticated models. Chuang, Kantor and Kardar (CKK) [16| 
have introduced a successful description of unbiased homopolymer translocation, assuming, 
essentially, that the polymer's evolution does not take it far from its equilibrium states. 
Astonishingly simple as it may sound, translocation is, then, ruled by the diffusion exponent 
of the polymer center of mass, as if there were no blocking membrane. 

We are interested to provide a general kinetic description of polymer translocation and, 
in particular, of the CKK scaling results. Elementary stochastic events of the underlying 
Markov chain processes are given by the translocation of individual monomers. As we show 
in this work, the CKK dependence of the translocation time with polymer size is indeed 
recovered by the Markov chain modeling, after finite-size corrections are properly eliminated 
- an issue not addressed by the original CKK approach. 

This paper is organized as follows. In sec. II, we write down the general Chapman- 
Kolmogorov equation which describes polymer translocation as a Markov chain, and discuss 
its exact asymptotic solution for the probability of complete translocation. The CKK picture 
is then taken into account as a way to devise expressions for the transition probabilities of 
individual monomer translocation events. The main point in our work is that the CKK 



scaling exponents may be considerably affected by finite size effects, which are partially 
encoded in corrections to the polymer mobility exponent. 

In sec. Ill, we compare results obtained from the numerical solution of the Chapman- 



Kolmogorov equation with the ones from previous detailed Langevin simulations 11 ]. 



As polymer translocation has been conjectured to be subdiffusive 17| and non-markovian 



18 



19| . we address, in sec. IV, a critical analysis of these issues, finding support, from 
further Langevin simulations, for the pertinence of the Markovian framework. In section V, 
we summarize our results and point out directions of further research. 

II. CHAPMAN-KOLMOGOROV APPROACH 

Our essential aim is to model polymer translocation as a discrete Markov stochastic 
process. Assume, for a homopolymer of length N, that at an arbitrary time instant there 
are n monomers on the trans side of the membrane and N — n monomers on the cis side, 
as depicted in Fig. 1. Let p n and q n be the probabilities that correspond to cis — > trans 
and trans — > cis monomer transitions, respectively. It is not necessary to have p n + q n — 1; 
actually, the polymer may get momentarily stuck around the membrane with probability 

1 - Pn ~ Qn- 

An object of central interest is the probability P(n, t) of finding n monomers on the trans 
side of the membrane at time t, once the initial probability vector P(n, 0) is given (t is an 
integer valued variable). In numerical experiments, it is usual to have P(n, 0) = 5(n,N/2), 
meaning that the initial state has half of the monomers on each side of membrane. The 
Chapman-Kolmogorov equation for the translocation process can be readily written as 

P(n, t + l) = q n+1 P(n + 1, t) + p n -iP(n - 1, t) 
+(l- Pn -q n )P(n,t) , (1) 

where we take < n < N, with p_ 1 = p = q = qw+i = Pn = <1n = 0, so that translocation 
proceeds until n = 0orn = N(n = and n = N are absorbing states of the stochastic 
process). 

The important task now is to find expressions for the transition probabilities p n and q n . 



In the CKK description of translocation 16( , it is assumed that the translocation time t^ of 



a polymer of size A^ is of the order of the Rouse relaxation time [20| , that is, the time spent 
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FIG. 1: The trans and cis sides of the membrane have, respectively, n and N — n monomers. 



by the center of mass of a free polymer to diffuse along its own gyration radius Rq ~ N u . 
Writing the difusion constant of the center of mass as D = D/N s , where D is the diffusion 
constant of a single monomer, we will have Rq ~ Dt, and, thus 



r N 



jy5+2v 



(2) 



The CKK expression tn ~ N 1+2l/ follows by the substitution of 5 = 1 (as it applies for free 
polymers in the absence of hydrodynamical couplings) in ([2]). In our discussion, however, we 
take 5 to be an adjustable parameter, which may slightly depart from unit for small chains 
due to pore/membrane interactions. An estimate of 5 can be obtained by the use of the 
polymer generalization of the usual Einstein relation which establishes the proportionality 
between the mobility and the diffusion coefficients in brownian motion 21(. We recall that 
while the mobility can be computed from the velocity response of the polymer to external 
forces, the diffusion coefficient refers to the polymer evolution in the absence of any external 
perturbations. In this respect, we point out that recent Langevin simulations 22] indicate 
that the mobility computed for driven translocation scales as N~ 5 , with 5 = 0.81 ± 0.04 for 
polymer sizes N < 256. 

It is important to note, before proceeding, that the original master equation approach 



introduced by Muthukumar 23[], which is based on a relaxation-to-equilibrium formalism, 
yields, very generally, the scaling law r ~ iV 2 for translocation. This result is a consequence 
from the fact that the transition probability ratio q n /p n would be given, according to the 



detailed-balance hypothesis of Ref. |23|, by 

"1 1 



— = exp<{ (7- 1) 

Pn 



n N — n 



(3) 



which is close to unit for 1 < n < N (above, 0.5 < 7 < 1 is the scaling exponent used in 
the definition of the partition function, Z n ~ n 7_1 , for a polymer chain anchored at one 



endpoint to a hard wall |23l . l24j]). The challenge we are faced, thus, is how to obtain the 
anomalous scaling predicted by the CKK picture from a markovian line of argumentation. 
A kinematical solution of this problem (i.e., not based on free-energy considerations) will be 
provided by an alternative choice of the transition probabilities p n and q n . As a bonus, we 
will be able to address the role of finite size effects in the translocation process. 

For a heuristic derivation of p n and q n , assume that the time lapses for single trans — > 
cis and cis — > trans monomer translocations are proportional to r n jn and r^_ n /{N — n), 
respectively. This implies that the transition probabilities satisfy the balance constraint, 

q n — = p n — = constant . (4) 

n N — n 

Substituting © in (gj), it follows that 

c 

Pn 



c 

,<5+2j/-1 ' 



Qn - n s+ 2v _i > ( 5 ) 



where < c < 1 is an arbitrary constant. As we will see in the next section, solutions of the 
Chapman- Kolmogorov equation (0Q) with transition probabilities (jSJ) will lead to transloca- 
tion times which are not given by (J2J), due to the existence of finite-size effects, a main point 
of attention in our work. We find that (|2J) is recovered only in the asymptotic limit of large 
polymer chains. 

The probability of complete translocation to the trans side of the membrane, P(N, 00), 
known as the "splitting probability" in the theory of stochastic process, is available in closed 
analytical form [25] . Since this quantity depends on the initial number of trans monomers, 
n, we write P(N, 00) = P(n), with 



Ml- 1 



i + Eg n; = i | 
i + eS 1 nj=i | 



We have performed Langevin simulations in order to check if the above exact splitting 
probabilities are in fact pertinent. The results, which support the viability of the Markov 
chain approach to polymer translocation, are reported in sec. IV. 

III. SCALING REGIMES 

The main difficulty in carrying out realistic Langevin numerical simulations of polymer 
translocation is usually related to limitations in the polymer sizes (typically, present com- 
puter desktop resources allow one to work finely with a few hundred monomers). It turns 
out that scaling results, like Eq. fl2]), are unavoidably affected by finite size scaling effects, so 
that even if statistical error bars are taken into account, predictions and observations may 
not satisfactorily match. 



We address in the following a comparison between Markov chain modeling resu 



ts and 



the ones obtained through the extensive Langevin simulations reported in Ref. 11] . The 
general strategy is to numerically solve the Chapman-Kolmogorov equation for a given set of 
transition probabilities p n , q n (corresponding to different translocation regimes) and several 
polymer sizes N, always taking the initial probability vector as P(n,Q) = 8(n,N/2). 

Let r be the translocation time for an individual realization of polymer translocation. It 
is clear that 

oo 

r = j}l - 6(a - r)] , (7) 

s=l 

where B(s — r) is the Heaviside step function of s — r (we are using the convention O(0) = 0), 
so that the mean translocation time can be written simply as 

oo 

r = (r)=X)[l-<e(*-r))] 

8=1 

OO 

= £[1 - P(N, a) - P(0, a)] , (8) 

8=1 

where the above averages are taken over the ensemble of polymer translocation realizations. 
In (jSJ), P(N, s) + P(0, a) is the probability, to be obtained from the numerical solution of the 
Chapman-Kolmogorov equation, for the occurrence of complete translocation up to time s; 



it is a monotonically increasing function of s, which approaches unit for s — )■ oo. For the sake 
of clarity, note that if F s is the probability that complete translocation takes place for the 
first time at time s, then P(N, s) + P(0, s) = YH>=o F 9 ' (recall that translocation is modeled 
here as a stochastic process with absorbing boundary conditions). Since "Y^ = qF s > = 1, we 
get 1 - P(N, s) - P(0, s) = J27=s+i F s'> and > therefore, from ©, 



oo oo 



r= EE F ^E sF - ( 9 ) 

s=0 s'=s+l s=l 

which is the more familiar (but unpractical, for our purposes) way of writing r as a mean 
first passage time. 

It suffices, for excellent numerical convergence, to retain in (JBJ) all the contributions which 
have P(N, s) + P(0, s) < 1 — 10~ 6 (results have precision better than 10~ 3 %). The scaling 
exponent a in r ~ N a is then determined by straightforward linear regression on log-log 
plots. 

We assume, in the quasi-equilibrium regimes discussed here, that both sides of the mem- 
brane are characterized by the same coil-globule transition temperature Tg, i.e., in a real 
experiment the solvent qualities would be the same on both sides of the diathermic mem- 
brane. There are, roughly, three possible translocation regimes depending on the equilibrium 
temperature T: (i) T < Tg, (ii) T = Tg and (in) T > Tg. We model these regimes, respec- 



tively, by gyration radius exponents 



21 



26 



27J (i) v = 1/3 (globule phase), (ii) v = 1/2 



(0-point), and (iii) v = 0.588 (self avoiding walk). 

We have found through numerical solutions of the Chapman-Kolmogorov equation that 
the translocation time scales with the polymer size in all of the above situations. In order 
to compare our results with the ones from Langevin simulations [ll| we have limited the 
number of monomer sizes up to N = 500. The scaling exponents and translocation time 
profiles for different translocation regimes are shown in Table I and Fig. 2, respectively. 
The finite size corrections are found to be in good agreement with the Langevin simulation 
results if one takes 5 ~ 0.88 ± 0.03. We note that the discrepancy for the T = Tg regime is 
probably due to lack of precision in locating the #-point in Langevin simulations (no error 



bars are reported in Ref. ll( for T = Tg). In passing, we call attention to the fact that as T 
drops much below Tg, it has been reported, in the same referred work, that the dynamical 
exponent ai (defined in Table I) stops decreasing and eventually grows to values larger than 
2. It is possible that such a non-monotonic behavior of q^ has to do with known deviations 



of Rouse's theory at low temperatures [281 ]. 
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FIG. 2: The translocation time r is found to scale with the polymer's size N as r ~ N a . We 
report scaling profiles for 50 < N < 500 and T < Tg (squares), a = 1.739(2); T = Tg (circles), 
a = 2.059(1); T > Tg trans (triangles), a = 2.231(1). We have considered 5 = 0.88. 



Temperatures 


OLth 


«L 


OiM 


T <Tg 


5/3 


1.74(3) 


1.75(3) 


T = Tg 


2 


2.2(?) 


2.07(3) 


T>Tg 


2.176 


2.23(3) 


2.24(3) 



TABLE I: Comparison between scaling exponents. Here, we assume otth = 1 + 2z/ |l6l ] for the 
quasi-equilibrium cases, while a^ and ctjf are determined through Langevin simulations [11] and 
the Markov chain approach (with 5 = 0.88 ± 0.03). Both ai and olm are evaluated for polymer 
sizes in the range 50 < iV < 300. The error bar for the Langevin simulation with T = Tg is not 
available. 



The crucial point here is that due to finite size effects, the scaling exponents evaluated 
from Langevin simulations could seem to contradict the theoretical CKK predictions, even if 



error bars are taken into account in numerical evaluations. The translocation time exponent 
for a polymer of size N can be written, in general, as 



a = 6(N) + 2v + f(6(N) + 2v, N) , 



(10) 



with 5(N) < 1, f(5(N) + 2v, N) > 0, and 



lim S(N) = 1 , 

AT— ►oo 

lim f(6(N) + 2i/,N) = . 

AT— >oo 



(11) 



We expect that 5(N) and f(S(N) + 2u, N) are relatively slow functions of N. 

For a fixed value of 5, we may compare evaluations of the scaling exponents obtained from 
the Markov chain approach for progressively larger polymers, with those proposed within 
the CKK framework. The finite size scaling results are shown in Fig. 3, where we take 5=1 
just for the sake of illustration. One finds, in fact, convergence towards the conjectured 
values of the scaling exponents. 
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FIG. 3: Finite size scaling behavior of a — (1 + 2v) (we have taken 5 = 1). Corrections to the 
scaling exponents a.th = 1 + 2v vanish as a power law N~P with scaling exponents /3 = 0.1418(5) 
(squares), j3 = 0.1348(4) (circles), and /3 = 0.1327(3) (triangles). 



Our data is consistent with a — otth ~ A r_/3 , where (3 ~ 0.13 — 0.14. This scaling be- 
havior implies that the translocation time can be written, taking into account subleading 
corrections, as 

T „ N 1+2u [l + cN~P log(iV)] . (12) 

The theoretical computation of (3 is an important challenge deserved for further studies. 
It is worth emphasizing that a complete discussion on the universality of the translocation 
exponents should necessarily come along with finite size scaling arguments. 

IV. EVIDENCE OF MARKOVIAN BEHAVIOR 

The Markov approach of this work stands in contrast to previous claims that memory 
effects are relevant in the phenomenon of polymer translocation. We intend here to go deeper 
in this issue, first through a critical account of the non-markovian line of thought, and then 
by pointing out further results which provide clear support for a markovian description of 
polymer translocation. 

Let s(t) be the number of monomers on a given side of the membrane at time t, in the 
course of unbiased polymer translocation. In the wake of the CKK's picture of translocation, 
it has been conjectured that s(t) should be modeled as a non-markovian process. The chain 
of ideas which would suggest the non-markovian character of polymer translocation can be 
summarized as follows: 

(i) Defining the time-dependent variance of the translocation coordinate s(t) as A(t) = 
(s(t) 2 ) — (s(t)) 2 , one expects that A(r) ~ iV 2 holds for a polymer of size N, where r is 
the translocation time. This result would naturally follow from the anomalous diffusion law 

2 

A(i) ~ tT+2l7; 

(ii) Since v > 0.5 in the polymer high-temperature phase, the anomalous diffusion expo- 
nent 2/(1 + 2v) < 1 indicates that polymer translocation is a non-markovian process; 

(iii) Monte Carlo simulations of polymer translocation, as performed within the frame- 
work of the Bond Fluctuation Model (BFM) [29|, are able to reproduce the anomalous 
diffusive profile advanced in (i). 

It is important to emphasize that there is a large room for questioning the relevance of 
the above statements, as we point out below, in respective order. 



10 



(i') There is no universal relation between the dynamic exponents associated to the 
translocation time r and the diffusion law. They are completely independent as a mat- 
ter of principle. Actually, we have checked from solutions of the Chapman- Kolmogorov 
equation ([T]) with the transition probabilities (J3J (which lead, as discussed in the previous 
sections, to the CKK translocation exponent 1 + 2v in the limit of large polymers) that most 
of the translocation process is described by the normal diffusion law, where A(£) is a linear 
function of time; 

(ii') Anomalous diffusion is by no means a sufficient condition for non-markovian behavior. 
There are abundant and important examples in the literature of either super or subdiffusive 



markovian processes [30|. We recall the theory of Levy flights, just to quote a celebrated 
instance of Markov processes which have anomalous diffusion exponents [311 ]: 

(iii') Even though Monte Carlo algorithms are not able in general to address dynamical 
aspects of statistical systems, the BFM has been accomplished as a very useful tool for the 



study of dynamical phenomena in polymer physics [32]. It is clear that the all the credit 
for the BFM's approach relies on the comparison of its predictions with the outcomes of 
real and numerical experiments. Despite the popularity of the BFM strategy in polymer 
physics, it is not obvious a priori if under the specific boundary conditions related to polymer 



translocation, BFM will provide physically meaningful results [33 1. 

Regarding (iii'), in particular, our attention is drawn to recent extensive Langevin sim- 
ulations of unbiased polymer translocation which have been performed in connection with 
the problem of anomalous diffusion 19j . Its is claimed, in that work, that on the basis of 
standard log-log plots, anomalous diffusion is verified with dynamic exponent 2/(1 + 2v) at 
intermediate time scales. However, we have found that if the very same data is represented 
in linear scales, as in Fig. 4, an excellent and unique linear fit for A(£) holds both for small 
and intermediate time scales. The Langevin simulations of Ref. [19J seem to give support, 
actually, to normal diffusion in the form A(£) = a + bt (we note that it is a straighforward 
exercise to derive arbitrary a and b coefficients from a simple collored gaussian noise version 
of the Ornstein-Uhlenbeck process, which is essentially markovian for time scales larger than 
the noise correlation time). 

We have also carried out Langevin simulations in order to investigate the splitting prob- 
abilities predicted by Eq. fl6]). While we have not produced a large number of realizations 
(compared to the 5000 of Ref. [19|) which would allow us to compute the diffusion exponent 

11 
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FIG. 4: The numerical data listed in Fig. 4a of Ref. [19[] is alternatively plotted in linear scales. 
The linear fit strongly suggests that for these time scales we have A(t) = a + bt, with a = 0.62 and 
6 = 2.3 x 1(T 5 . 

with reasonable precision, we have collected 560 complete translocation processes, which 
are enough to address a numerical test of Eq. ([6]). Our polymer has N = 50 monomers 
which translocate through a pore denned at the center of the membrane, taken as an 80 x 80 
monoatomic lattice. The pore is created by the remotion of a single atom of the membrane. 

nn 

Following the usual Langevin modeling prescriptions [9|, 111], the mononomer- monomer or 
monomer-membrane interactions are described by the Lenard- Jones potential 

4e[(a/r) 12 - (arf] + e , r < 2 1 / 6 cx 
, r > 2 1 / f V . 

Consecutive monomers in the polymer chain interact, additionally, through the finitely ex- 
tensible nonlinear elastic (FENE) potential, 

1 



Ulj{t) 



(13) 



U F (r) = --kR 2 \n[l - (r/i? ) 2 ] . 
Monomers evolve, then, according to the Langevin equations, 
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FIG. 5: The empirical (black squares) and analytical (red circles) probabilities (both denoted here 
by P(n)) of cis — > trans complete translocations are compared. The variable n stands for the initial 
number of monomers in the trans-side of the membrane. The analytical probabilities are computed 
from (ED and (H with 5 = 0.88 and v = 0.588. 



where r,- 



r\,|, £ is the dissipative constant and Fi(t) is the gaussian stochastic force 



which acts on the monomer with label i, 



([n ■ Fi(t)][h' ■ F 3 (t')}} = 2n ■ tik B Tti5 l0 5{t - if) . 



(16) 



Above, n and h' are arbitrary unit vectors, and k B and T are the Boltzmann constant and 
the temperature, respectively. By means of a suitable regularization of the stochastic force, 
we have implemented a fourth-order Runge-Kutta scheme for the numerical simulation of 
ffT6|) . Our simulation parameters are: e = 1.0, o = 1.0 (a is also identified to the membrane 
lattice parameter), £ = 0.7, k = 7e/a 2 , R = 2a, k B T = 1.2e. The simulation time step is 
taken as 3 x 10 -3 £lj, where tu = \Jma 2 /e is the usual Lenard- Jones time scale. 

The initial configuration of the polymer has n monomers on the trans-side of the mem- 
brane. Translocation is allowed to occur only after thermal equilibrium is reached for the cis 
and trans sectors of the polymer. The empirical probability of translocation to the cis-side of 

13 



the membrane is then measured for n = 5, 10, 15, ...,45. For each n, we consider an ensemble 
of 70 complete translocations. 

In order to compare Eq. (JHJ) with the empirical probability computed from the Langevin 
simulations, it is necessary to realize that the Markov approximation is likely to hold for 
the consecutive translocation of small monomer clusters, rather than for the consecutive 
translocation of individual monomers. We expect, thus, that even though there may be 
strong correlation effects between the translocation of consecutive monomers, these correla- 
tions become small from cluster to cluster (it is reasonable to assume that the cluster size 
is proportional to the polymer persistence length). We have found, as it is shown in Fig. 
5, that a suggestive comparison between the empirical and analytical probabilities can be 
obtained considering clusters with the size of 5 monomers. More precisely, for a given value 
of n (the initial number of trans-monomers) a red circle is plotted in Fig. 5 with coordinates 
(n, P(n/5)) where P(n/5) is the splitting probability evaluated for a polymer which contains 
iV/5 = 10 monomers. 

The empirical error bars in Fig. 5 were evaluated by the partition of each translocation 
sample (which has 70 elements) into 7 sub-sets. The theoretical error bars, on the other hand, 



follow from elementary statistical considerations and are given by \J P(n)(l — P(n))/70. 

We note that due to slow variation of the translocation dynamic exponents with the 
polymer size, the existence of the small monomer "markovian clusters" , as introduced above, 
has a neglegible effect on the results established in Sec. III. 

V. CONCLUSIONS 

We have discussed in this paper a general kinetic model of polymer translocation, within 
the framework of markovian stochastic process. In order to comply with the standard picture 
of translocation [16j, we have put forward a Chapman-Kolomogorov equation with transition 
probabilities p n and q n , which turn out to depend as power laws on the number of cis (or 
trans) monomers in the polymer chain. We have also established a closed analytical expres- 
sion for the probability of complete polymer translocation in terms of arbitrary p' n s and q' n s. 
The Chapman- Kolmogorov equation is then numerically solved for quasi-equilibrium regimes 
where the CKK picture of translocation is assumed to work. We have been able to find good 



agreement with the scaling results derived from previous realistic Langevin simulations 11], 



14 



only if a non-trivial scaling relation, due to finite size effects, is taken into account for the 
diffusion constant of the polymer center of mass, a phenomenon observed in simulations by 
Bhattacharya et al. |22|]. Our results indicate, thus, that finite size effects, which have been 
so far left to a secondary role in most of the polymer translocation literature, should have 
fundamental importance in the analysis of real or numerical translocation experiments. It 
would be interesting, in particular, to revisit Langevin simulations of polymer translocation 
under various solvent qualities, having in mind the relevance of finite size corrections. 

Subleading corrections to the predicted asymptotic scaling profiles are found to have a 
slow decay, a fact that could explain some of the controversy on the precise values of the 
scaling exponents, and the issue whether they are actually universal or not. The reasonable 
precision associated with not very large polymers in Langevin simulations (see the third 
row in Table I) could suggest at first that universality or the CKK picture needs revision. 
However, we regard this apparent difficulty as a peculiar effect due to subleading scaling ex- 
ponents and to corrections on the polymer mobility exponent. We point out, in this respect, 
that the amplitude of subleading corrections can be sensitive to specific modeling details, an 
additional complication factor in the empirical evaluation of universal translocation scaling 
exponents. 

The Markov chain modeling strategy can be used, in principle, as a valuable tool in 
further real or numerical translocation experiments. Once a set of transition probabilities 
p n and q n is determined in an experimental study of translocation, a comparison between 
the markovian expression ((6]) and the observed frequency of complete translocation as a 
function of the initial number of trans monomers can be carried out, in order to reveal (or 
not) the existence of memory effects in the translocation process. Actually, we have been 
able to perform a successfull test of Eq. ([6]) with the help of Langevin simulations and 



have noted that a re-interpretation of recent Langevin simulation results [19| suggests that 
polymer translocation diffuses in a normal way (as a linear function of time), which is in 
complete agreement with the Markov chain picture addressed here. 

This work has been partially supported by CNPq and FAPERJ. 
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